*=c 


COMPUTATION  OF  THE  TURBULENT  MIXING 
IN  CURVED  EJECTORS 


Change  of  Address 


Organizations  receiving  reports  on  the  initial  distribution  list  should 
confirm  correct  address.  This  list  is  located  at  the  end  of  the  report. 
Any  change  of  address  or  distribution  should  be  conveyed  to  the  Office 
of  Naval  Research,  Code  211,  Arlington,  VA  22217. 

Disposition 

When  this  report  is  no  longer  needed,  it  may  be  transmitted  to  other 
organizations.  Do  not  return  it  to  the  originator  or  the  monitoring 
office. 

Disclaimer 


The  findings  and  conclusions  contained  in  this  report  are  not  to  be 
construed  as  an  official  Department  of  Defense  or  Military  Department 
position  unless  so  designated  by  other  official  documents. 

Reproduction 

Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 


w-  0Ni^|:R212-249-2F  J/ _ foi 1 

4.  TlTLg  (**d  Subiitt*)  _  __ 

^  '  j  COMPUTATION  OF  THE  JURBULENT ^MIXING  INCURVED  j 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  I'nr.l  fWtum  t>*l*  hn  frill) 

Tiffi  REPORT  DOCUMENTATION  PAGE  befomk^compuktino^form 

■  <  W  t.t  QM  h«  r_»tCIRlCMT‘S  CATALOG  NUMAEA 


0NR\CR212-249-2F 


I  A.  TITLE  (ltd  Suiitttt) 


»/  TYPE  of  REPORT  >  RERIO 

f  Final  Jfepwrfc. 

(  1  Sept  J»78  -  31  A 

«■  RCRFORMIHGlOAO.  REPORT  Nuua 

NR8^H-ll7 

TTTOtTaC?  OR  GRANT  KUM 


EJECTORS. 


J  E.  F.'Schum  \  ^  - -  -i 

\  P.  H.l  Bev  ilaquaj  //<-  V  Nl>jDQ14-78-C-p557f<< 

j  S.  V./Patankar  J _ \Z^  L  ~ - ' 

».  PfAfOAM'NC  ORGANISATION  NAME  ANO  ADORES*  10.  PROGRAM  ELEMENT.  PROJ 

.......  .  „  AREA  A  WORK  UNIT  MUMOI 

/^Rockwell  Ir.ternationali^vijjuiat'icn  %  ^ _ _  ^^im  s-  / 

"Rorth  American  Aircraft  Divtsieis/  //v  T  RF4  1&7m3j61  J ^-r' 

Co’.umbi’S ,  0N»  4321o 

M.  CONTROLLING  OFFICE  NAME  AND  AOORESS  ^ -  5.  allSt**1*  T"*-'**  C~~  "  ' 

Oft  ice  of  Javal  Research  fit'  Apr®  ®80 

Vehicle  Technology  Program  (Code  211)  (  J  TTinJGIF'Wmror- 

Arlington,  VA  22217  N'— 35 

U."'MArt7i  Jning  agency  mane  a  aqones i/t/ ftw  Ca*tr*Uuis  otUc»)  »».  security  class.  (»i  uu* 

iJt/  j  Unclassified 

- - J  IS«.  OCCLASSiriCLTidN'/'boWNGNAdtNO" 

^  SCHEDULE 

la”  WSTHlljklTlON  VfATeMeNT  f»l  UU.  "  ■  — — 


Approved  for  public  release,  distribution  unlimited. 


IT.  USTNteuTlON  STATEMENT  f«T  ia»  #Numi  mHh^ to  RterA  iA.  Itaihwi  a—  EmhU 


hi-  Su^lEUEnTANy  notes 


It  At  »  ■OHStS  teaa<iHM»  s*  w.«M  mO*  Uwmuai 

ejectors 

V/STOL  aerodynasies 
turbulence 

viseid/iaviseid  flov  interactions 


An  AU>* 

aircraft  propulsion 
thrust  augmentation 
jet  flap  diffuser 


A  numerical  method  has  been  developed  to  calculate  the  effeets  of  curvature  00 
the  turbulent  mixing  within  a  thrust  augmenting  ejector.  Curvature  eodifiea- 
e ion  are  included  in  both  the  kinetic  energy  and  dissipation  equations  of  the 
turbulence  model.  The  elliptic  effects  of  curvature  on  the  aean  flow  are 
determined  by  combining  an  elliptic  solution  for  the  preasure  field  with  a 
conventional  marching  solution  for  the  velocity  field.  Predictions  of  the 
velocity  profiles  for  jets  and  curved  duct  flows  are  compared  with  experi¬ 
mental  data.  Computation  of  the  curvature  effeets  are  shown  for  an  ejector 


00  (*anMTS  1473  COITION  0*  1  NOVAS  iSOOSOlETC 


_ _ UNCLASSIFIED _ 

security  ctnuncaTioa  o*  this  rage  j» 


c/2 


SCCumfry  Cl  ASSiriCATlOH  or  THIS  MA0KtH7.*n  nmtm  Fnffd) 


UNCLASSIFIED 


having  Coanda  jets  and  one  central  Jet* 

T 


a  5tfu»i!>  CW  Q»  TaUi  P*6{*4(«a  tin 


TABLE  OF  CONTENTS 


Sec  tion  Title 


P£££ 


TABLE  OF  CONTENTS . ill 

LIST  OF  FIGURES .  iv 

LIST  OF  TABLES .  v 

NOMENCLATURE  .  vi 

INTRODUCTION  .  1 

GOVERNING  DIFFERENTIAL  EQUATIONS  .  4 

Mean  Flow  Equations . . . .  4 

Turbulence  Model.  .....................  5 

DIFFERENCE  EQUATIONS  AND  SOLUTION  PROCEDURES  .  9 

General  Iteration  Procedure  ................  9 

Grid  and  Difference  Equations  ...............  10 

Pressure  Scheme  ......................  lb 

Grid  Generation  .  ..............  18 

RESULTS  AND  DISCUSSION  .  20 

Flow  in  a  Straight  Duct  ..................  20 

Flow  in  a  Curved  Duet  ...................  22 

Flow  in  an  Ejector . .  2b 

CONCLUSIONS .  30 

REFERENCES .  31 


LIST  OF  FIGURES 


Figure  No.  Title  Page 

1  XFV-12A  Integrated  Lift/Propulsion  System.  .  .  1 

2  General  Solution  Procedure  .  9 

3  Grid,  Velocity  and  Control  Volume  Notation  .  .  11 

4  Slab  Continuity. . . . .  15 

5  Orthogonal  Squares  Used  in  Ejector  Flow  Field 

Analysis  ...  .  .......  19 

6  Comparison  of  Predicted  and  Measured  Velocity 
Profiles  for  Fully  Developed  Turbulent  Flow  in 

a  Channel . 20 

1  Comparison  of  Predicted  and  Measured  Velocity 

Profiles  for  a  Plane  Turbulent  Jet  in  a  Co- 
f  lowing  Stream . 21 

8  Comparison  of  Predicted  and  Measured  Peak  Vel¬ 

ocities  in  a  Turbulent  Wall  Jet  in  a  Co-flowing 
Stream  ....................  22 


9  Comparison  of  Predicted  Velocity  Profiles  with 

Mathematically  Exact  Theoretical  Solutions  for 
Fully  Developed,  Laminar  Flow  in  2-Dimensional 
Ducts  with  and  without  Curvature  .......  23 

10  Comparison  of  Predicted  and  Measured  Velocity 
Profiles  for  Turbulent  Flow  in  a  Rectangular* 

Curved  Duct.  .................  24 

11  Comparison  of  Predicted  and  Measured  Pressures 

for  Turbulent  Flow  in  Rectangular*  Curved  Duet  25 

12  Comparison  of  the  Cross  Stream  U  Velocity  Pro¬ 
files  Obtained  from  the  Parabolic  Program  and 

the  Curved  Ejector  Program  ..........  2b 

13  Variation  of  the  Cross  Stream  Pressures 

Obtained  from  the  Curved  Ejeetor  Program  ...  V 


tv 


NOMENCLATURE 


a 

A 

Area 

b 

BB 

C 


cp  c2,  c n 

d 

D 

E 

C 

k 

l 

io 

L 

N 

P 


»i»  *o 

He 

S 

u  V 
u* 

♦  * 

v* 

V 

§P 

VOL 


U 


convection  and  diffusion  term  defined  by  Equation  (18) 
convection  and  diffusion  term  defined  by  Equation  (18) 
area  used  in  the  convection  and/or  diffusion  flow 
width  of  jet  nozzle 

residual  flow  in  the  continuity  equation,  Equation  (32) 
convection  rate.  Equation  (17) 
turbulence  constants 

diffusion  rate  defined  by  Equation  (16) 
diffusion  or  substantial  derivative 
nodal  point  E  (Figure  3) 

rate  of  generation  of  turbulent  kinetic  energy,  Equation  (9) 
turbulent  kinetic  energy 

mixing  length  with  curvature  effects  included 

mixing  length  in  the  absence  of  curvature  effects 

width  of  a  channel 

nodal  point  N  (Figure  3) 

pressure  or  nodal  point,  P  (Figure,  3) 

distance  in  radial  direction  or  radius 

inner  or  outer  radius 

Reynolds  number 

distance  in  the  r®  direction  or  nodal  point  S  (Figure  3) 

Reynolds  stress  term 

turbulent  fluctuation  velocity  in  the  U  direction 
velocity  in  the  $  direction 
velocity  on  eenterline 

terms  used  in  the  tl  momentum.  Equations  (24)  and  (23) 
turbulent  fluctuating  velocity  in  the  V  direction 
velocity  in  the  radial  direction 

terms  used  in  the  V  momentum,  Equations  (27)  and  (28) 
volume 

turbulent  fluctuating  velocity  in  a  perpendicular  direction 
from  U  and  V 

nodal  point  U  (Figure  3) 


i 


X 

y 

a 

6 

( 

9 

k 

P 

*k 

r 

0 

W 


distance  from  nozzle  exit 

distance  measured  from  a  wall  (or  centerline)  in  a  channel 
(eg,  Figure  6) 

empirical  constant  used  in  Equation  (2) 

distance  used  in  diffusion  equation,  (Equation  16)  or  shear 
layer  thickness  (Equation  2) 

rate  of  dissipation  of  turbulent  kinetic  energy 

angle  used  in  curvilinear  coordinate  system 

bulk  viscosity 

turbulent  viscosity 

density 

Prandtl  number  for  use  in  the  turbulent  kinetic  energy 
equation  (Equation  10) 

Prandtl  number  for  use  in  the  turbulent  dissipation  equation 
(Equation  11) 

symbolic  notation  used  in  Table  1 
symbolic  notation  used  in  Table  1 
angularity  velocity 


Subscripts 

E 

in 

jet 

aax 

N 

out 

P 

S 

see 

U 


East 

in. 

je* 

maxiaTO 

North 

out 

at  eenter 
South 
secondary 
West 


superscripts 

u  for  U  aoaentua  equation 

v  for  V  aoeentue  equation 


vii 


■it  v  ii :  am  i  rriiii 


v.  wif'ik  Mi 


3 


INTRODUCTION 


The  U.  S,  Navy  is  considering  several  categories  of  V/STOL  aircraft 
for  use  with  smaller  carriers,  as  a  more  economical  means  of  maintaining 
sea  control.  One  way  of  obtaining  the  additional  thrust  required  by 
these  aircraft  is  to  divert  the  engine  exhaust  flow  through  a  thrust 
augmenting  ejector.  Such  an  ejector  is  a  kind  of  jet  pump  in  which 
entrainment  by  the  exhaust  stream  accelerates  a  larger  mass  of  air  drawn 
from  the  atmosphere.  By  Newton's  law  of  action  and  reaction,  the  ejector 
experiences  a  force  which  is  equal  but  opposite  to  the  momentum  change  of 
the  entrained  air.  A  more  complete  description  of  this  process  has  been 
given  by  Bevilaqua.l 


Combining  the  ejector  with  the  wing,  as  shown  in  Figure  1,  produces  an 
especially  effective  lift/propulsion  system.  A  V/STOL  aircraft  having  an 


Figure  1.  XFV-UA  integrated  i-irt/Propulsioo  System 
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ejector  wing  converts  smoothly  from  jet-borne  hover  to  wing-borne  flight, 
because  the  wing  lift  is  increased  during  conversion  by  the  jet  flap 
effect  of  the  ejector  exhaust  stream.  In  addition,  separate  reaction 
jets  are  not  required  for  control  during  hover;  control  forces  may  be 
generated  by  differential  action  of  the  ejectors.  To  demonstrate  these 
advantages,  the  North  American  Aircraft  Division  oe  Pockwell  International 
is  developing  the  XKV-12A,  a  V/STGl  aircraft  which  has  ejectors  in  the 
wing  and  canard. 

Although  analytic  methods  for  predicting  ejector  wing  performance  are 
necessary  for  conceptual  design  studies  and  to  reduce  test  requirements, 
a  theory  for  predicting  both  the  lift  and  thrust  of  ejector  wings  has  not 
been  developed.  There  are  methods  for  calculating  the  increments  of 
aerodynamic  lift  and  pitching  moment  induced  by  an  ejector  of  given 
thrust.  These  are  based  on  the  now  classical  vortex  sheet  model  of  the 
pure  jet  flap  devised  by  Spence. *  linearized,  thin  airfoil  theories  were 
developed  for  the  ejector  wing  by  Chan*  and  Woolard4,  who  added  a  sink 
on  the  upper  surface  of  a  jet-flapped  wing  to  represent  entrainment  by 
tlie  ejector.  Wilson-*  extended  this  approach  by  analysing  complete 
geoatetries ,  including  thickness  and  camber  as  well  as  deflection  of  the 
jet  wake.  Mere  recently,  Mendenhall  and  McMillan6  have  begun  investiga¬ 
tion  of  three  dimensional  effects. 

These  methods  have  been  useful  in  performing  parametric  analysis  and 
for  cctaputing  load  distributions,  such  as  surface  pressures  and  hinge 
moments.  However,  neither  the  thrust  augmentation  ratio  nor  the  initial 
thrust  angle  is  predicted,  in  current  practice,  experimental  data  is 
used  to  specify  the  variation  in  these  parameters  during  conversion  fr^m 
hover  to  conventional  flight.  Such  an  empirical  approach  is  only  useful 
for  small  variations  from  the  original  data  base,  so  that  there  is  no 
procedure  for  evaluating  significant  design  changes  or  new  configurations. 

For  static  ejectors,  hevilaqua  and  Gedoode*  developed  a  method  of 
predicting  the  thrust  augmentation  by  iterating  between  a  viscous  solution 
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for  he  turbulent  mixing  within  the  ejector  and  an  inviscid  solution  for 
the  flow  outside  the  ejector.  The  purpose  of  this  report  is  to  describe 
a  viscous  solution  developed  for  use  in  a  similar  analysis  of  the  ejector 
wing.  Calculating  the  turbulent  mixing  is  more  difficult  in  this  case, 
because  curvature  effects  become  significant  as  the  ejector  is  deflected 
aft.  In  the  next  section,  the  effect  of  curvature  is  discussed  and 
the  governing  equations  are  presented.  The  numerical  scheme  devised  to 
solve  these  equations  is  described  in  the  following  section.  In  the  last 
section,  the  calculated  results  for  some  simple  flow  geometries  are  com* 
pared  with  available  data,  and  a  sample  calculation  for  a  typical  ejector 
configuration  is  shown. 


GOVERNING  DIFFERENTIAL  EQUATIONS 


Mean  Flow  Equations 

The  performance  of  the  ejector  wing  depends  primarily  on  the  thrust 
augmentation  ratio  and  the  thrust  angle,  however,  computing  these  para¬ 
meters  is  made  difficult  by  the  complexity  of  the  ejector  flow  field. 

There  are  interacting  regions  of  turbulent  and  irrotational  flow,  subject 
to  lateral  straining  and  streamwise  curvature.  The  curvature  produces 
strong  lateral  pressure  gradients,  which  distort  the  inlet  flow  and 
deflect  the  thrust  vector  from  the  mean  diffuser  angle. 

These  elliptic  effects  due  to  curvature  cannot  be  calculated  with 
classical  parabolic  solution  procedures, in  which  upstream  influences 
are  neglected  and  the  equations  are  solved  by  marching  through  the  ejec¬ 
tor  in  the  streamwise  direction.  While  a  fully  elliptic  solution  would 
include  upstream  influences,  the  requirements  for  increased  computer 
storage  and  time  make  such  an  approach  prohibitively  expensive.  Instead, 
we  have  adapted  a  method  devised  by  Pratap  and  Spalding9* *0  for  economi¬ 
cally  calculating  upstream  influences  in  duct  flows.  Since  there  is  a 
primary  direction  of  flow  (through  the  duct)  the  elliptic  effects  are 
primarily  due  to  transmission  upstream  by  the  pressure  field.  This  type 
of  flow  falls  between  parabolic  flows,  in  which  there  are  no  upstream 
influences,  and  elliptic  flows,  in  which  pressure,  convection,  and  dif¬ 
fusion  transmit  influences  in  every  direction.  From  a  computational 
viewpoint,  duct  flows  of  this  type  may  be  classed  as  partly  elliptic  or 
partially  parabolic. 

The  equations  that  govern  the  flow  in  the  ejector  are  derived  from  Reynolds' 
equations  for  turbulent  flows,  by  neglecting  streamwise  diffusion  and  upstream 
convection.  Thus,  the  effect  of  a  local  disturbance  is  transmitted  upstream 
by  presure  only;  downstream  by  pressure  and  convection;  and  across  the 
flow  by  pressure,  convaction,  and  diffusion.  In  cylindrical  coordinates 
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(r,  9)  the  governing  equations  take  the  form 


Conservation  of  Mass 


au  ev  ,  V 

rdQ  ar  r 


=  0 


Streamwise  Momentum 


n  +  v  M  ■  uv  I  SP_  du'v'  2u  V 
rae  Sr  r  r  "  p  rde  "  ar  -  r 


Radial  Momentum 


(1) 


av  av  u2  i  ap 

rdG  ar  r  p  d r 

These  equations  do  not  describe  the  effect  of  curvature  on  the  flow. 
The  "extra"  terms  originate  in  the  choice  of  coordinate  system  and  do  not 
imply  any  additional  physical  processes.  The  Coriolis  acceleration, 

UV/r  =  o»V  and  the  centrifugal  acceleration  U*/r  =  u>U,  arise  from  rotation 
of  the  velocity  vector  wj respect  to  the  coordinate  system.  Similarly, 
the  extra  stress  term  is  a  consequence  of  having  chosen  a  cylindrical 
volume  element;  that  is,  since  the  sides  of  the  element  are  not  parallel, 
the  shear  stress  has  a  component,  u'v'/r,  in  the  streamwise  direction. 
Thus,  the  use  of  cylindrical  coordinates  is  a  convenience  in 
this  case,  and  the  additional  terms  merely  account  for  the  difference 
between  the  actual  accelerations  of  the  fluid  and  those  which  appear  to 
take  place  in  the  curvilinear  coordinate  system. 

Turbulence  Model 

The  effect  of  curvature  is  to  increase  the  turbulent  mixing  rate.  The 
first  attempts  to  model  this  effect  were  straightforward  extensions  of 
mixing  length  theory,  Prandtl**  and  Townsend**  developed  relations  to 
increase  the  mixing  length  as  the  radius  of  curvature  decreased.  These 
had  the  general  form 

I  *  /o  U  +  *  “)  (2) 
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in  which  l  is  the  mixing  length,  b  is  the  shear  layer  thickness,  and  ct 
is  an  empirical  constant.  This  approach  significantly  underpredicts 
the  change  in  the  mixing  rate.  Bradshaw^  improved  this  model  by  making 
the  mixing  length  increase  in  proportion  to  the  strain  rate  due  to  curva¬ 
ture, 

i  =  l0  ^1  +  a  Qu/57)  •  (3) 


More  recent  proposals  have  been  based  ou  the  turbulence  kinetic  energy 
equations.  According  to  the  usual  eddy  viscosity  assumption,  the  turbu¬ 
lent  stress  is  first  expressed  in  terms  of  a  turbulent  viscosity,  /l^,  and 
the  mean  strain  rates 


-  = 


(4) 


Again,  no  new  physical  process  is  implied  by  the  extra  strain  rate, 

U/r.  Following  Launder  and  Spalding1^*  the  turbulent  viscosity  is  assumed 
to  depend  on  two  parameters:  the  turbulence  kinetic  energy,  k,  and  its 
rate  of  dissipation,  «.  From  dimensional  arguments,  the  expression  for 
the  turbulent  viscosity  is 


**t 


(5) 


in  which  c/^  is  an  empirical  constant. 
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Under  the  partially  parabolic  flow  assumptions,  the  curvature  equations 
for  the  three  components  of  the  turbulence  kinetic  energy  take  the  following 
forms : 


Change  «  Production  -  Dissipation  +  Diffusion  +  Rotation 
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In  contrast  to  the  mean  flow  equations,  these  turbulence  equations 
contain  additional  terms  which  do  represent  a  physical  process  caused  by 
the  curvature.  Since  the  extra  terms  have  opposite  signs  in  the  equations 
fjr  u'  and  v',  the  effect  of  curvature  may  be  interpreted  as  an  eddy  rota¬ 
tion  which  transfers  energy  from  the  streamwise  component  of  the  turbulence 
to  the  transverse  component.  This  interpretation  of  the  curvature  terms 
was  suggested  by  Hunt  and  Joubert.^ 

Summation  of  these  three  equations  gives  the  equation  for  the  turbulence 
kinetic  energy,  k  =  ylu'^  4-  v'2  +•  w'^)  used  by  Launder  and  Spalding;^ 

',uS  +  <,vi-D  +  G  •  *<  <7> 


Because  the  rotation  simply  transfers  energy  without  loss  between  the 
turbulence  components,  ''-irvature  is  seen  to  have  no  effect  on  the  total 
turbulence  energy.  Thus,  the  classical  turbulence  energy  equation  cannot 
be  used  to  determine  curvatur~  effects. 


Instead,  we  have  used  a  hybrid  equation,  which  simulates  the  actual 
physical  mechanism,  .he  turbulence  kinetic  energy  is  used  to  provide 
an  estimate  of  the  average  intensity  of  the  velocity  fluctuations  which 
cause  the  turbulent  mixing.  However,  in  thin  shear  layers  it  is  the  v* 
component  of  the  turbulence  whleh  produces  most  of  the  mixing.  There¬ 
fore,  an  additional  production  tern  was  added  to  the  turbulence  energy 
equation  to  represent  the  increase  of  v*  due  to  curvature.  The  equation 
for  the  turbulence  kinetic  energy  therefore  becomes 


pu 


0  k 
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(8) 


in  which  C  is  the  production  of  turbulence  kine-k  energy  which  occurs 
in  the  streamwise  component.  It  has  the  fora 


The  expression  for  the  diffusion  of  energy  by  the  turbulent  fluctuations 

is 
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(10) 


An  equation  derived  by  Chambers  and  Wilcox^  was  used  to  compute 
the  dissipation  of  the  turbulence  energy.  It  contains  the  same  curvature 
convection  as  the  turbulence  energy  equations 


pu  +pv  a?  =  (ciG  "  c2 p°  k 


dt  \Ot  dx  ) 


This  turbulence  model  has  five  empirical  constants.  According  to  the 
recommendations  of  Launder  and  Spalding^  the  following  values  of  the 
constants  were  used  in  this  analysis. 


c^  C|  C2  Ok  Of 

0.09  1.44  1.92  1.0  1.3 

Kor  boundary  layer  flows.  Chambers  and  Wilcox^7  suggest  using  9/2, 
rather  than  the  factor  of  2,  in  the  curvature  term.  The  effect  of  chang¬ 
ing  this  value  was  not  found  to  be  significant. 


DIFFERENCE  EQUATIONS  AND  SOLUTION  PROCEDURES 
General  Iteration  Procedure 

The  equations  involved  are  a  coupled  set  of  nonlinear,  partial  dif¬ 
ferential  equations.  These  equations  are  transformed  into  a  corresponding 
set  of  finite  difference  equations  that  necessitate  an  iterative  type  of 
solution  to  obtain  the  velocity  and  pressure  field.  The  iterative  approach 
is  needed  because  the  coefficients  in  the  basic  differential  equations  also 
involve  the  velocities  and  pressures  being  sought.  Figure  2  briefly  depicts 
the  overall  iterative  approach  and  this  is  followed  by  a  more  detailed  dis¬ 
cussion  of  the  method  of  solution. 


Step  l.  Initialization 

Estimate  entire  pressure  and  velocity  field  at  problem  start. 

Step  2.  Matching  Solution  far  Velocities 

Solve  u- momentum  equations  for  downstream  U  velocities. 

Correct  U's  to  satisfy  overall  mass  continuity. 

Solve  continuity  equation  for  V  velocities. 

Solve  turbulence  equations  for  P  s«d  t  yielding  Mi. 

Repeat  for  next  downstream  slab. 

Step  1.  Solution  for  Entire  Pressure  Field 

Solve  Poisson's  equation  for  pressures  at  eaeh  node. 

Step  4.  Convergence  Check  of  Solution 

Return  to  Step  2  if  pressures  and  velocities  are  not  converged. 


Figure  2.  General  Solution  Procedure 


The  major  improvement  in  this  approach  is  the  manner  in  which  the  pres¬ 
sures  are  calculated.  In  the  earlier  method  of  Bevilaqua  and  DeJoode,^ 
a  single  pass  or  march  through  the  ejector  was  made.  Pressures  were 

stored  in  a  single  array  in  the  program  and  downstream  pressure  effects 

9  10 

could  not  be  transmitted  upstream.  In  the  work  of  Pratap  and  Spalding  * 
many  passes  were  made  through  the  duct  and  on  each  pass  the  downstream 
pressure  perturbation  was  transmitted  upstream  over  an  incremental  step 
distance.  Convergence  to  the  final  pressure  field  was  relatively  slow. 

This  is  not  the  case  for  the  program  used  herein.  Pressures  are  stored 
in  a  two-dimensional  array.  The  equations  used  to  calculate  all  of  the 
nodal  pressures  are  solved  in  a  simultaneous  manner  ("sweeping").  Any 
downstream  effect  is  transmitted  all  the  way  upstream  on  each  pass  through 
the  duct.  Convergence  is  much  faster. 

Grid  and  Difference  Equations 

The  notation  for  the  velocities,  U  and  V,  and  their  control  volumes  fol¬ 
low  the  same  general  procedure  used  by  Pratap,  et  al*°.  The  control  volumes 
are  staggered  as  shown  in  Figure  3  and  the  subscripts  P,  N,  S  and  E  refer 
to  the  particular  variables  directed  toward  these  nodal  points.  The  nota¬ 
tion  N,  S,  E  and  W  designate  the  points  to  the  North,  South,  East  and  West. 

Using  this  notation  the  continuity  and  momentum  equations  can  be 
expressed  in  a  simplified  form. 

Mass  Conservation: 

(pUp)  •  Area  -  (pUg)  ♦  Area  +  (pVp)  •  Area  -  (pV^)  •  Area  «  0 
U  Momentum: 

APUP  s  +  AvUv  +  As^s  +  *  Area  +  Source  *  Volume  (12) 

V  Momentum: 

ApV p  *  A^V,,  +  a\v  +  AvV  ♦  (P  -  p  )  •  Area  +  Sooree  •  Volume 
rr  cv  v  v  5  \r 
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(a)  Nodal  Point  and 

(b)  Control  Volume  for 

Velocity  Notation 

Continuity,  Pressure 

lr 

k  and  € 

isa 


(c)  Control  Volume  for  U  Velocity 


kviffi 


(d)  Control  Volume  for 
V  Velocity 


Figure  3,  Grid,  Velocity  and  Control  Volume  Notation 

where  the  coefficients,  A,  represent  the  combined  effects  of  convection 
and  diffusion  while  the  source  terms  include  the  effects  of  curvature,  e 
The  area  terms  correspond  with  the  incremental  volume  surface  over  which 
the  velocity  or  pressure  is  aeting.  To  greatly  simplify  the  computer 
program  logic,  Pa  tanka  r  ^  Iras  suggested  that  both  the  momentum  and 
turbulence  equations  be  expressed  in  the  simplified  fora: 


puH  +  ^vlf 3  ^  (ril)+ Source  (13) 

I 

\ 

Convection  Diffusion  j 

[ 

i 

in  which  the  symbol  0  can  represent  U,  V,  k  or  €  as  shown  in  Table  I,  j 


Table  I.  Generalized  Equation  Parameters 


Equation 

l 

r 

Source 

U  Momentum 

U 

d P 
e)S 

U  d  v  puy 

•  ~i  Tx  -  V 

V  Momentum 

V 

0 

e>P 

dr 

r 

Turbulent  Kinetic  Energy 

k 

Mt/ak 

G  - 

Rate  of  Dissipation  of 
Turbulent  Kinetic  Energy 

£ 

/a€ 

(c,c  -  CjP«  -  f  »t  ^ 

If 

Using  the  source  terra  notation  of  Table  I,  Equation  (15)  can  be 
expressed  in  a  general  finite  difference  forra  as: 


Ap  0jj  +  Ay  0w  +  As  0s  +  Souree  *  Volume  (14) 

in  which  the  subscripts  refer  to  the  nodal  points  in  Figure  3.  To  insure 
mathematical  stability,  the  “Source"  terns  in  the  above  equation  ahould 
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contain  only  the  positive  values  of  those  shown  in  Table  I  while  the 
negative  source  terms  should  be  substituted  into  the  Ap  relation! 


t 


i 


< 


Ap  =  Ajg  +  Ay  Ag  -  Source  •  Volume  (15) 

The  terns  A  in  this  equation  represent  a  combination  of  diffusion,  d, 
and  convection,  C.  For  example,  the  diffusion  from  point  P  to  N 
(Figure  3)  for  use  in  the  U  momentum  equation  can  be  described  as: 


a 


Area 


&N  5P 


(16) 


in  which  the  interfacial  area  term,  Area,  is  located  at  a  distance,  bjj, 
from  the  N  grid  point  and  distances  bp,  from  the  P  grid  point.  Although 
the  turbulent  viscosity  is  shown  in  Equation  (16),  laminar  values  can  be 
used. 

There  is  also  a  convective  flow,  from  nodal  point,  P,  to  N. 

»  P •  Area  *  V  (1?) 

To  insure  a  more  rapid  and  stable  convergence  in  |ae  iterative  solution, 
Pataakar*7**8  has  developed  the  following  relations  for  use  in  Equations 
(14)  and  (15). 


******  [*  °] 

where 

a*  »  dN  (l  -  .1  ICjjl/djj)*  for  |c^  |  <  10  (IS) 

a*j  »  for  (Cji  10  djj 


\ 


.[]. 


in  which  the  brackets,  L  J,  signify  the  largest  of  the  two  enclosed  fac¬ 
tors.  In  like  manner,  the  combined  diffusion  and  convection  term,  Ag, 


between  nodes  S  and  P  are; 


AS  =  as  +  [cs,  oj 


The  coefficient  Aw  in  Equations  (14)  and  (15)  represents  the  convection 
on  the  west  side  of  the  control  volume; 

Ay  =  P  *  Area  •  Uw  (20) 

In  practice,  the  values  of  the  diffusion  and  convection  terms  are  averaged 
over  two  adjacent  control  volumes.  It  is  recotmoended  that  reference  to 
Patankar ^ ^  be  made  for  a  more  detailed  description. 

The  marching  technique  (Figure  2)  starts  with  the  first  slab  and  proceeds 
downstream,  slab  by  slab.  For  example,  the  Ap,  A(j,  Ag,  Aw  along  with  the 
source  terras  are  calculated  for  each  cell  in  the  particular  slab  considered. 
The  unknown  U  velocities  (Equation  (14)  with  0  set  to  U)  in  the  slab  are 
calculated  using  a  tridiagonal  matrix  algorithm  with  the  velocity  at  the 
wall  set  to  zero  for  the  wall  friction  case.  These  U  velocities  generally 
do  not  satisfy  the  requirements  of  continuity  over  the  entire  slab  (unless 
the  convergence  of  all  the  momentum  and  turbulence  equations  is  obtained) 
and  an  adaptation  of  the  method  of  Sparrow,  et  al^  is  used  to  “adjust** 
the  0  velocities  to  satisfy  slab  overall  continuity. 

The  V  velocities  are  then  calculated  for  each  cell  in  the  slab.  Figure  4, 
starting  with  the  lowest  cell,  using  the  basie  continuity  equation  (con¬ 
stant  density). 


VGOt  * 


Ufa  «  At*ajn  -  ti0iu  .  AceaQut  *  V|B  AreatR 
A?eaoui 


in  which  the  indicated  areas  ?orrespond  with  the  velocities.  It  wiil  be 
shown  later  that  the  converged  values  of  V  depend  on  both  continuity  and 
the  radial  momentum  equation  (Equation  l). 
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Figure  4.  Slab  Coot inuity 

rose  the  particular  slab,  the  values  of  U  and  V  are  continually  being  up¬ 
dated.  Using  the  "asst  reeent“  values,  the  turbulent  kinetic  energy,  k, 
and  dissipation,  «,  equations  are  solved  (Table  i  and  Equations  8  through 
ill  in  a  similar  eanne?  as  the  U  astsentwa  equations,  the  newly  calculated 
values  of  k  and  €  foe-  eae».  ceil  are  then  substituted  into  Equation  (5)  to 
detemiae  the  turbulent  viscosity,  jjt,  for  the  particular  cell. 

This  entire  procedure  is  repeated  for  each  succeeding  slab  until  the  aug- 
esenter  exit  is  reached.  At  this  point  in  the  iteration  process,  an  estimate 
of  the  entire  flow  field  of  u  and  V  velocities  is  known,  although  these 
values  eay  not  be  the  final,  converged  values. 


Pressure  Scheme 


Before  describing  the  detailed  equations,  it  should  be  pointed  out 
that  there  are  two  distinct  improvements  in  the  determination  of  the 
overall  pressure  field  when  compared  with  that  done  earlier  (eg,  Bevilaqua 
and  DeJoode^).  The  first  is  that  the  pressure  field  in  the  program 
described  in  this  paper  involves  an  elliptic  solution  in  that  down¬ 
stream  pressures  are  transmitted  to  upstream  regions.  The  second 
major  improvement  is  the  detachment  of  the  pressure  equation  from 
the  momentum  equation  in  the  solution.  In  the  earlier  approach 
a  pressure  correction  was  developed  and  was  used  to  perform  two  tasks-- 
correct  the  pressure  and  correct  the  velocity.  This  resulted  in  improved 
velocity  distributions  but  relatively  bad  pressure  fields  in  that  con¬ 
vergence  was  slow.  In  the  new  improved  approach,  two  separate  equations 
are  used,  one  to  obtain  the  pressure  (solution  to  Poisson  type  equation) 
and  the  other  to  correct  the  velocity.  Test  eases  have  shown  that  despite 
the  extra  numerical  computations  needed  in  the  new  procedure,  there  is  a 
jv  tu  iu  percent  saving  In  computer  time. 


The  following  paragraphs  describe  the  procedure  used  for  the  calcula¬ 
tion  of  the  pressure  field.  These  pressures  are  needed  for  substitution 
back  into  the  momentum  and  turbulence  equations  to  obtain  an  updated  vel¬ 
ocity  field  ( figure  1). 


The  d  momentum  equation  has  the  form 

ApUp  *  *  AyOy  ♦  AgU^  ♦  Source  *  Vo  l  e  Area  *  (Py  »  Pp)  (22) 


in  which  the  area  is  the  inter facial  area  perpendicular  to  the  line  extend¬ 
ing  from  node  point  Py  to  Pp  (Figure  5).  Ap  is  defined  by  Equation  (1$) 
and  includes  the  convect iua/dif fusion  terms  A$«  Ay,  and  Ag.  ke-arranqing 
Equation  (22)  leads  to; 


Up 


*[ 


Ajjtfjj  ♦  Aytly  ♦ 


A^tfg 


♦  Source 


.  Veil 

- T 


Area 


-  Pp) 


(25) 


lu 


(24) 
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and  replacing  the  bracketed  term  by  Up  results  ini 


Up  «  Up  +  ^I£a  (py  .  pp) 


In  like  manner  (with  different  values  of  area) 

ue  "  "e  +^r  “V  •  pe> 


The  V  momentum  equation: 


aNvN  +  AWVW  +  ASVS  +  Source  *  Vol]  Area 

AP  J  *  ~  {Ps  '  Pp) 


can  also  be  simplified  to: 


(25) 


(26) 


VP  55  %  +  (?S  “  (2?) 

and 

vn  a  vN  «►  Aga  (pp  -  pji)  (28) 

The  Vp  and  VN  relations  represent  the  radial  Boeentusa  equations.  With  the 
substitution  of  Equation  (23)  into  the  continuity  equation  for  the  pressure 
control  volume  (Figure  i  ): 

(pU  Arca)y  -  (pU  Area)g  *  (pV  Area)8  -  (pV  Area);;  »  0  (29) 

there  is  obtained 

ApPp  «  A^PS  ♦  AyPy  +  A§P$  ♦  AgPg  *55  (jy) 

in  which 

A?  *  A$  ♦  Ay  ♦  A§  ♦  Ag  (51) 

and 

65  3  (P  Are.  dp)  .  ( p  Area  %)  ♦  (pArea  Vp)  -  (  pArea  0j»)  (52) 


I 
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With  convergence  the  value  of  BB  tor  each  cell  over  one  entire  field 
approaches  zero. 

There  are  five  unknowns  in  Equation  (32);  namely,  the  pressures 
at  the  center  of  the  cell  and  the  pressures  of  the  four  adjacent 
cells.  When  the  number  of  cells  in  the  flow  field  is  large,  it 
is  expedient  to  solve  for  the  pressures  using  a  block  relaxation  technique 
followed  by  the  use  of  the  tridiagonal  matrix  algorithm  sweeping  from  left 
to  right,  then  right  to  left  followed  by  sweeps  in  the  vertical  directions. 
For  this  study  Patankar^^  developed  a  "block  relaxation"  technique  that 
was  shown  to  greatly  accelerate  convergence.  In  brief,  a  constant  pres¬ 
sure  correction  is  added  to  each  slab  prior  to  the  sweeping  action  type 
computations. 

The  newly  calculated  pressures  are  then  substituted  back  into  the 
momentum  and  turbulence  equa '  .ons  (Step  2  in  Figure  2)  and  the  entire 
procedure  is  repeated.  Convergence  is  obtained  when  further  calculations 
show  that  the  values  of  .elocity,  k,  €,  and  pressure  do  not  change  within 
chosen  tolerances. 

Grid  Generation 

In  older  to  calculate  the  velocity  field  in  a  non-uniform  area  duct, 
as  in  an  augmenter,  it  is  necessary  to  first  generate  a  nearly 
orthogonal  grid  system,  forming  a  set  of  "curvilinear  squares"  to  which 
the  finite  difference  equations  can  be  applied.  For  accuracy,  a  large 
number  of  such  squares  or  cells  is  desired.  In  the  absence  of  a  compre¬ 
hensive  computer  program  to  perform  this  operation,  the  flow  field  in  an 
augmenter  was  first  subdivided  into  a  coarse  grid  of  squares  by  a  trial- 
and-error,  hand  drawing  technique  (eg,  Schneider‘S).  Richardson^ 
indicated  agreement  between  graphical  and  exact  analytical  solutions  of 
tne  order  of  1.0  percent  for  simple,  classical  problems.  The  coarse  grid 
for  the  augmenter  analysis  is  presented  in  Figure  5.  About  midway 
through  the  augmenter,  the  "squares"  were  extended  into  rectangles  for 
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Figure  5.  Orthogonal  Squares  Used  in  Ejector  Flow  Field  Analysis 


core  enience  in  the  computations.  A  computer  program  was  written  to 
further  subdivided  the  coarse  grid  into  a  finer  grid,  consisting  of  a 
total  of  3375  celis,  45  in  the  cross  stream  direction  and  75  cells  in 
the  flow  direction.  The  necessary  geometric  information  for  each  of 
these  cells  was  recorded  on  computer  output  tape  for  read- in  to  the 
viscous  program.  Results  of  the  augmenter  flow  field  calculations  as 
well  as  some  simpler  test  cases  are  presented  in  the  following  section 


RESULTS  AND  DISCUSSION 


The  comparison  of  experimental  data  with  analytically  predicted  values 
from  a  computer  program  can  be  used  to  assess  the  accuracy  of  the  mathe¬ 
matical  modeling  of  the  complex  flow  phenomenon.  In  this  section  comparisons 
arenadefor  flow  in  straight  ducts,  as  well  as  in  curved  ducts  in  order 
to  demonstrate  the  cross  stream  (elliptic)  pressure  effect.  Flow  in  an 
augmenter  is  also  studied. 

Flow  in  a  Straight  Duct 


A  comparison  of  predicted  and  measured  velocities  for  fully  developed 
turbulent  flow  in  a  2-dimensional,  parallel  plate  channel  is  shown  in 
Figure  6.  Good  agreement  with  a  mean  line  through  the  data  of  Laufer^ 
is  obtained. 


A  similar  comparison  is  shown  in  Figure  7  for  the  case  of  a  plane 
turbulent  jet  in  a  co-flowing  stream,  the  data  being  that  of  Weinstein, 
et  al.  The  initial  jet  velocity  is  twice  that  of  the  secondary  flow,  a 
value  which  is  not  too  far  different  from  that  for  augmenters. 


Figure  7.  Comparison  of  Predicted  and  Measured  Velocity 

Profiles  for  a  Plane  Turbulent  Jet  in  a  Co- flowing 

Stream 


At  distances  to  60  jet  widths,  (the  maximum  investigated  in  the  tests), 
the  agreement  is  good.  The  slight  differences  at  the  outer  edge  of  the 
jet  might  be  attributed  to  the  use  of  a  "thick  walled"  nozzle  in  the 
experiment. 
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Figure  8.  Comparison  of  Predicted  and  Measured  Peak 
Velocities  in  a  Turbulent  Wall  Jet  In  a  Co¬ 
flowing  Stream 

The  predicted  deeay  of  the  maximum  velocity  of  a  wall  jet  in  a  co- 
flowing  stream  is  compared  in  Figure  8  with  the  measured  results  of 
Cartshore  and  Newman^  *  Good  agreement  is  also  shown. 

Flow  in  a  Curved  Duet 


Although  the  program  was  primarily  written  to  study  turbulent  mixing 
in  ducts  (augmenters)  the  program  also  has  the  capability  for  calculating 
laminar  flow  In  ducts  with  and  without  curvature.  The  effect  of  curvature 
on  fully  developedi  laminar  flow  Is  shown  in  Figure  9  for  a  duct  with  a 
high  degree  of  curvature. 
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Figure  Comparison  of  Predicted  Velocity  Profiles  with 
Matheaat leal iy  Exact  Theoretical  Solutions  for 
Fully  Developed,  laminar  Flow  in  2-Dimensional 
Duets  with  and  without  Curvature 

Computer  predicted  values  are  io  good  agreement  with  the  exaet  theoreti¬ 
cal  solution,  developed  by  So-5  .  As  a  matter  of  interest,  the  theoreti¬ 
cal  exaet  solution  is  also  shown  for  the  ease  where  there  is  no  eurvature. 
For  laminar  flew  it  can  he  seen  that  ihe  e fleet  of  eurvature  is  small, 
witch  is  not  the  ease  for  turbulent  flow. 

The  turbulent  flow  data  of  Eskinaci  and  Teh26  are  used  here  for 
comparison  with  predieted  values.  The  test  configuration,  Figure  10,  con¬ 
sists  of  4.aS  meters  (lb  ft)  straight  duct  preceeding  the  eurved  duct  which 
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of  the  straight  and  curved  sections,  predicted  velocities  agrse  well  with 
a  mean  curve  through  the  data.  At  the  170°  station  and  inner  surface 
(y/L  <  .3)  predicted  velocities  are  larger  than  measured  values  while  at 
the  outer  surface  (y/L  >  .7)  the  predicted  values  are  smaller.  The  maxi¬ 
mum  difference  in  the  calculated  and  measured  velocities  is  of  the  order 
of  10  percent.  Calculated  pressures  across  the  duct  are  in  excellent 
agreement  with  measured  values  (Figure  11).  The  effect  of  curvature  on 
the  velocity  profile  is  seen  to  be  much  larger  for  turbulent  flow  than 
laminar  flow,  when  he  results  of  Figures  9  and  10  are  compared. 
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Figure  11.  Coaparison  of  Predicted  and  Measured  Pressures  for 
Turbulent  Flow  in  Rectangular.  Curved  Duct  , 


The  ejector  model  (Figures  5  and  12)  has  a  throat  width  of  0.122  meters 
(0.4  ft)  and  a  distance  or  length  from  the  throat  to  the  exit  plane  of 
0.268  meters  (0.88  ft).  A  characteristic  geometric  parameters  is  the  ratio 
of  this  length  to  the  ejector  throat  width.  For  this  configuration  the 
ratio  of  length-to-width  is  2.2,  The  diffuser  half  angle  is  10°.  The 
center  jet  consists  of  a  plane  2-dimensional  slot  jet  with  a  total  width 


Figure  12.  Comparison  of  the  Cross  Stream  U  Velocity  Profiles 
Obtained  from  the  Parabolic  Program  and  the  Curved 
Ejector  Program 
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of  0.762  cm  (0.3  in)  while  the  Coanda  gap  is  .191  cm  (0.075  in).  Nozzle 
exit  velocities  are  315.2  meters/sec  (1034  ft/sec).  The  ratio  of  the 
radius  of  the  Coanda  surface  to  the  Coanda  gap  is  approximately  17.6  at 
the  jet  nozzle.  Calculated  velocities  and  cross  stream  pressure  dis¬ 
tributions,  obtained  from  the  curved  ejector  program,  are  shown  for  3 
stations.  At  the  throat.  Station  A,  the  secondary  velocity  increases  from 
85.3  meters/sec  (280  ft/sec)  near  the  center  jet  at  y  ■  1.37  cm  (0.045  ft) 
to  94.5  meters/sec  (310  ft/sec)  near  the  Coanda  jet  at  y  ■  5.18  cm  (0.17  ft), 
the  increase  due  to  curvature.  Over  the  same  region  the  pressure  in  the 
secondary  flow  decreases  from  -4310.  to  -7180.  pascals  (-90  to 


Pressure 

(Pascals) 


figure  13.  Variation  of  the  Cross  Stream  Pressures  Obtained 
from  the  Curved  Ejector  Program 
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-150  psfg)  as  shown  in  Figure  13  causing  a  V  component  or  flow  toward 
the  Co and a  surface  and  hence  an  increase  in  the  U  velocities. 

At  station  B,  the  secondary  velocities  decrease  with  distance  from  the 
centerline,  in  contrast  to  that  at  station  A.  Curvature  effects  at  sta¬ 
tion  B  are  small  (large  radius  of  curvature)  as  indicated  by  the  negligi¬ 
ble  cross  stream  pressure  gradient  in  Figure  13  at  this  station.  At  the 
exit,  the  two  jets  are  nearly  merged  and  cross  stream  pressure  gradients 
are  likewise  negligible. 

For  comparison  purposes,  the  calculated  velocities  obtained  with  the 
parabolic  flow  program  (Beviiaqua7)  are  also  presented  in  Figure  12  in 
which  the  secondary  mass  flow  rate  is  about  4  percent  larger.  This  small 
difference  in  flow  rates  should  have  little  effect  on  the  overall  trends. 

In  the  parabolic  program,  curvature  in  the  flow  field  for  ejector  con¬ 
figurations  eannot  be  accounted  for  nor  can  the  downstream  pressure  effects 
have  any  influence  on  the  upstream  pressures  and  velocities  (parabolic 
type  solution).  As  a  result,  calculated  secondary  velocities  are  constant 
at  any  axial  station,  as  shown.  However,  in  the  eurved  ejector  programs 
the  individual  nodal  pressures  and  corresponding  velocities  are  inter¬ 
related  over  the  entire  flaw  field.  In  the  comparison  of  eenter  jet 
velocity  profiles,  it  should  he  noted  that  the  maximum  or  centerline 
velocities  are  larger  for  the  paraboiie  program  and  this  can  be  attributed 
to  the  relative  axial  location  of  the  jet  nogsie.  In  the  parabolic  program, 
the  location  of  the  eenter  jet  must  be  in  the  same  plane  as  the  Coanda 
noeeles.  In  the  curved  ejector  program,  the  center  jet  is  located  upstream 
of  the  Coanda  noazle  plane  as  indicated  in  Figure  Su  With  the  center  jet 
no?ale  farther  upstream,  the  downstream  jet  profile  will  be  wider  and 
have  a  .ewer  peak  velocity  than  that  shown. 

In  the  comparisons  shown  in  Figure  12,  the  Coanda  jet  ooasles  are  in 
the  same  axial  location  in  the  ejector.  Kesuite  fur  the  curved  ejector 
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program  show  a  faster  decay  of  the  maximum  wall  jet  velocity  with  a 
corresponding  thickening  of  the  wall  jet.  This  is  due  to  the  addition 
of  curvature  effects  because  at  station  A  there  is  a  pressure  decrease 
of  approximately  4790  Pascals  (100  psf)  through  the  wall  jet  (from  y  •  5.2 
to  5.94  cm).  The  effect  of  curvature  on  Coanda  entrainment  should  be 
studied  in  more  detail,  in  particular  for  Coanda  surfaces  having  a  small 
radius  of  curvatures.  Although  a  symmetric  ejector  was  analyzed  in  this 
paper,  to  permit  comparison  with  the  results  of  parabolic  ejector  pro¬ 
grams,  asymmetric  configurations  may  also  be  analyzed  with  the  scheme 
described  here. 
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CONCLUSIONS 


A  viscous  turbulent  fluid  flow  program  has  been  developed  to  extend 
the  knowledge  of  the  complex  flow  phenomena  occurring  in  ejectors  during 
the  conversion  or  transition  portion  of  the  flight.  For  these  conditions 
curvature  of  the  flow  through  the  ejector  can  cause  a  large  pressure 
gradient  across  the  duct,  normal  to  the  main  flow  direction.  Because  of 
this  and  its  unknown  effect  on  thrust,  the  curvature  relations  have  been 
included  in  both  the  mean  flow  and  turbulence  equations  to  assess  the 
effect  of  curvature  on  the  flow  field. 

To  establish  a  confidence  level  in  the  program,  comparison  eases  were 
studied  in  which  the  program  was  used  to  calculate  the  velocities  and 
pressures  in  straight  as  well  as  curved  channels,  and  for  wall  jets  and 
center  jets  located  in  a  eo-r loving  stream.  All  of  these  flow  eases  can 
occur  in  an  ejeetor.  Sueh  comparisons  showed  good  agreement  of  predicted 
and  measured  data.  For  a  static  ejector  the  pressure  gradients  across 
the  ejector  channel  were  shown  to  be  large.  Curvature  effects  in  the 
flow  resulted  in  an  increase  in  the  spreading  rate  (increased  entrainment) 
of  the  Coaada  jet. 

based  an  the  favorable  results  obtained  with  this  program,  it  should 
prove  very  useful  in  the  development  of  high  performance  augmenters,  in 
particular  when  this  viscous  flew  program  is  “inter -connected1*  with  an 
inviseid  program  used  far  the  calculation  of  the  velocities  outside  of 
the  ejector. 
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